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^ Abstract. We analyze a system consisting of an oscillator 

i ■ coupled to a field. With the field traced out as an environ- 

bX), ment, the oscillator loses coherence on a very short decoherence 

j> \ timescale; but, on a much longer relaxation timescale, predictably 

evolves into a unique, pure (ground) state. This example of re- 
coherence has interesting implications both for the interpretation 
of quantum theory and for the loss of information during black 
hole evaporation. We examine these implications by investigating 
the intermediate and final states of the quantum field, treated as 
an open system coupled to an unobserved oscillator. 
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1 Introduction 

1.1 Overview of decoherence 

Recent years have witnessed a significant increase of interest in the process of 
decoherence^, 0, §, [| - - the loss of quantum coherence suffered by a quantum 
system in contact with an environment. An environment consists of degrees 
of freedom which are coupled with, but not regarded as an integral part of the 
system. External variables are an obvious example, but even internal degrees 
of freedom may constitute an environment, if they cannot be followed by the 
observer. Many quantum systems are therefore subject to decoherence, and 
the phenomenon is thus widespread and important. 

For example, it has been demonstrated that collective observables of 
macroscopic quantum systems will lose quantum coherence very quickly by 
this means. This loss of coherence will proceed at very different rates, de- 
pending on the initial state of the system. Indeed, in the simple models of 
quantum apparatus proposed to describe the process of measurement, one 
can select an interaction Hamiltonian which commutes with the observable 
of the recording apparatus 0. The varying susceptibility of initial states to 
decoherence then allows one to model apparent collapse of the wave packet. 
The results can be taken to imply that different outcomes of a measurement 
are all present, but --in the language of Everett |J - belong to different 
branches of the universal wave vector. Their simultaneous detection is im- 
possible: decoherence leads to environment-induced supers election rules^, 
which effectively exclude a majority of states from the Hilbert space of the 
open system. In the context of the Many Worlds Interpretation (MWI) of 
quantum mechanics, environment-induced superselection supplies a preferred 
basis, which selects the "branches' into which the universal state vector is 
"splitting" . Decoherence can thus be thought of as a "missing link" between 
the quantum universe and classical reality, in that it provides the criterion 
for selection of preferred observables (such as position) while supplying an 
effective definition for classicality, as well as the rationale for the apparent 
"collapse of the wave packet" . 

The process of decoherence has also been studied in the somewhat less 
idealized, but still exactly solvable, model of an oscillator system coupled 
to a quantum field representing the environment. The evolution of this sys- 
tem is known as quantum Brownian motion^, |7], H], and it also exhibits 



environment-induced superselection [y]. One can describe superselection in 
the case of quantum brownian motion by appealing to the predictability 
siei>e(T(J - a formal implementation of the idea that the preferred quan- 
tum states will be stable (i.e., will minimize entropy production) in spite 
of the coupling to the environment. For example, the predictability sieve 
selects coherent states as the preferred states of an underdamped harmonic 



oscillator [|ITJ . Moreover, because decoherence occurs on timescales which 
are typically much smaller than a system's dynamical timescales (or even 
timescales associated with monitoring by an observer), one can convincingly 
argue that similar "nice" states will be singled out by environment-induced 
superselection in more realistic (and more complicated) situations. 

Much of the discussion of the implications of the decoherence process is 
based on the tacit assumption that "decoherence is forever" ; a concern is of- 
ten voiced that any sign of re-coherence would be trouble^]. In this paper 
we explore the problem of re-coherence, by exhibiting a model — quantum 
Brownian motion with a zero-temperature environment - - in which deco- 
herence happens quickly, but then gets "undone" slowly. Since decoherence 
and information are intimately related in quantum theory, it will turn out 
that this system can serve as an instructive toy model for the information 
problem in black hole evaporation. 

1.2 Outline of the calculation 

A brief and simplified preview of our calculation is in order. We will start our 
oscillator in a "Schrodinger Cat" state, a superposition of two well-separated 
coherent states. The initial state of the oscillator-field system will therefore 
beQ 

|* i ) = (c + |V + ) + a-|V-»|0>/wd- (1) 

(More complicated pure initial conditions can also be considered.) Interac- 
tion with the environment will (approximately, and for a few decoherence 
timescales) lead to 

\%) = {c+\if> + ) + c-\il;-.))\0)f idd ^ 

c + W + {t))\Q+(t))f ieUl + c-\il>4t))\$4t))f ield = |*(t)> , (2) 



^dn fact, we will consider states that are only approximately direct products of field 
and oscillator states; but the simplified outline presented in this Section still describes the 
essential physics involved. 



where 

($ + (t)|$_(t))«l. (3) 

Hence the density matrix of the oscillator will be given by 

p = |c+l 2 \Mt))(Mt)\ + M 2 IiM*)XiM*)I + 0«<MOI*-(*)» • (4) 

Thus, at this stage one might feel justified in "declaring victory": the de- 
coherence has happened, as the form of the density matrix of the system 
demonstrates. 

However, our system is a damped harmonic oscillator. In contact with 
the vacuum of a quantum field, it will slowly (on the relaxation timescale) 
approach the unique ground state regardless of the initial state. Hence, after 
a sufficiently long time, \^{t)) will approach 

|tt(oo)> =c + |0)|$ + ) + c_|0)|$_) 

= |0>(c + |* + >+c_|*_», (5) 

where |0) is — in our example — the ground state of the harmonic oscillator. 
Thus, decoherence seems to have "gone away": it did not prevent the state 
of the oscillator from re-cohering into a unique, pure state. Moreover, as a 
consequence of this recoherence the environment (field) has been put into a 
very awkward, pure "Schrodinger Cat" state of its own! 

Before this process has been completed, each of the two systems involved 
(the oscillator or the field, with the other traced out as unobserved) was in 
a mixed state, by virtue of the correlations between them. However, in the 
end these correlations have all disappeared. Or, to put it more precisely, 
the oscillator-field correlations have been "used up" to force the field into a 
highly non-trivial state. This state has the property (as one can anticipate 
intuitively, and as we shall prove in more detail below) that, in spite of its 
undeniable purity, it appears to be mixed when explored by approximately 
local measurements (limited in time and space to less than the duration of 
the recoherence episode). Thus, the information which appeared to be "lost" 
to observers who could access only one of the two systems (namely, the field) 
eventually re-emerges, but in a very obscure and hard-to-exhibit form. 

1.3 Analogy with black hole evaporation 

This sequence of events — possible rapid loss of information, and its eventual 
re-emergence after a long time, but in a form difficult to decipher - - is 



analogous to another interesting and fundamental process which has received 



a lot of attention in recent years: black hole evaporation |[13|, |14] , |15 
There, gravitational collapse rapidly increases the entropy of the Universe by 
the (usually large) difference between the entropy of the collapsing material 



and the final entropy of the black holeJT^], which is given (in bits) by the 
area of its horizon measured in units of square Planck length. 

The collapse seems to increase the entropy of the Universe, because the 
inside of the black hole horizon is inaccessible to external observers. Further- 
more, black hole evaporation puts the field outside the horizon into a mixed 
state which, when analyzed layer by layer, appears to contain approximately 
black body radiation, with entropy at least as large as the entropy shed by 
the black hole. There appears to be no information in the emitted radiation. 
So, when at the end of the process the black hole is gone (as seems likely), the 
entropy of the Universe is larger by at least the entropy increase which oc- 
curred during the collapse. This process can be analyzed in some detail in the 



case of the Witten black hole in 1 + 1 dimensional spacetime[|l5|, [16|]. Results 



in this case are usually cited (in spite of the facts that the analogy with the 
3 + 1 dimensional case is only partial, and that the calculation cannot really 
be carried out for the time when the black hole remnant ultimately disap- 
pears) as evidence for the black hole information paradox^H^ . This paradox 
is that the fundamental equations of gravitation and field theory seem to im- 
ply an irreversible increase of entropy in the process of collapse, despite the 
fact that they are themselves reversible. Reversibility cannot, of course, be 
established in this case simply by appealing to the dynamics, since Einstein's 
theory necessarily predicts a singularity inside the horizon, where the known 
laws of physics, including general relativity, are thought to fail. It is easy to 
imagine that the reversibility of the the whole process is an early victim in 



this failure [IS 



The decoherence cum re-coherence process described and analyzed in this 
paper supports an alternative view. One can imagine that the information 
lost beyond the black hole horizon eventually re-emerges, but not in any ob- 
vious form such that it could be detected by looking at "natural" observables 
(anything reasonably local, or at least confined to finite shells of the radiation 
emitted by the black hole). Rather, the information is re-emitted in a horri- 
bly "scrambled up" manner, where the state left behind after the black hole 
evaporates can still be pure (or at least no more mixed than the pre-collapse 
state), but this purity can only be revealed by measurement of some uncom- 



promisingly global observable, which coherently and simultaneously samples 
all of the emitted quanta. 

Motivated by this analogy, we shall exhibit such global observables, which 
can be computed exactly in our system by virtue of its linearity. It should of 
course be emphasized that this linearity that allows our model to be exactly 
solvable also makes it a rather distant analogue for the black hole evapo- 
ration process. It is precisely the inherent nonlinearity of general relativity 
which is responsible for the central singularity, the event horizon, the "no 
hair" theorems, and therefore for the unique value of the black hole entropy. 
Nevertheless, the complexity of the pure "global states" of the field gener- 
ated in our simple example suggests how the information can be preserved 
but remain "hidden" , and thus, suggests a possible resolution of the black 
hole information paradox. 

2 The calculation 

2.1 The model 

Our model consists of a simple harmonic oscillator Q coupled to a massless 
scalar field <p(x) in one dimension. To ensure that the energy is bounded 
below, we choose the following Lagrangian 0, with an ultraviolet cut-off: 

M n . 1 r°° 

L = -i?(Q 2 ~ tf Q 2 ) + - dx (d t <f>) 2 - (d^) 2 

Z Z J — oo 

/oo 
dxF(x)d t <f>{x) . (6) 

-oo 

The coupling constant g has dimensions such that we can define from it the 
frequency 70 = -f^-. Once renormalized, this frequency will correspond to a 
relaxation timescale. F(x) is a smearing function which implements the cut- 
off on the field-oscillator interaction. For our later convenience, we choose 
the particular form 
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dk — , cos kx 
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11 + k 2 

-K (T x) . (7) 



Kq(Tqx) is the modified Bessel function of order zero, which is concentrated 
within x < Tq 1 , and is a delta- function representation in the limit r — > oo. 
We quantize our model by defining the Hamiltonian operator 



H = - / rfx ([7r(x) + (?F(x)g] 2 + [^0(x)] 

Z J— oo v 



+^ (m^P 2 + Mofi 2 Q 2 ) . (8) 

Here we set % = c = 1, and introduce the convention that operators have 
circumflex accents, while c-numbers do not. P and jt(x) are the canonical 
momenta. 

In fact, this model is unitarily equivalent to several other much-studied 
systems, including even the free massless field. In writing (H), therefore, 
we are really choosing how the Hilbert space is to be divided into field and 
oscillator degrees of freedom. Our criteria for doing this in the way that leads 
to (H) are that we demand that an identifiable oscillator exists, that it be 
coupled locally (save for some UV smearing) to the field at the origin, and 
that the expected energies for direct product states of field and oscillator are 
finite. (The latter stipulation is needed in order for weak coupling to imply 
small entanglement between oscillator and field at low energies; it amounts 
to a demand for an ultraviolet cut-off.) 

To see this unitary equivalence, and to proceed in our calculation, we 
diagonalize (^j by defining the following normal modes, for u > 0: 

1 f<x> 

<p(x) = —= / du A w Uu{x) + B^sinux 

V7T JO 

tt(x) = —= / dull Vu{x) + Il5 sincux 

V7T JO 

i r°° 

Q = duH^q(u) 

y/7T JO 

1 f°° 

1 duj Ajp{w) ■ (9) 

7T JO 

The mode functions u^ and v U) and the co-efficients q(u) and p(co), 
are found by solving ordinary differential equations derived by iterating the 
Heisenberg equations of motion. One obtains 
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_u^_ 2 7o rw 2 
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4 7o r C(a;) 
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F(x) 
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-M u 2 q(uj) 



(10) 



We will soon see that u u {x) can be given in a much more transparent form. 
The important normalization co-efficient C(u) is defined so that 



C 2 (u) 



1 



u> 



n 2 



+ 



2 7o I> 2 \ 

nm+uj 2 )) 
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' 2 7o rga; 

v « 2 (rg + a^). 



2 + r 2 ] [cu 2 - (o + i 7 ) 2 ] [cj 2 - (n - i 7 ) 2 ] 
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(11) 



The convenience of the cut-off scheme defined in ([!]) lies in the fact that 
C 2 [lS) has only six simple poles, occurring in pairs of equal magnitude. This 
will make it easy to perform analytically several integrals that appear in the 
time evolutions discussed in later subsections. In Eq.(6) the new quantities 
r, fi,7 are modified versions of the frequencies r ,fi , 7 o, renormalized by 
the interaction. The renormalized frequencies may be expressed in terms of 
the bare ones only through a cumbersome (although analytically solvable) 
cubic equation. It turns out, however, that if we fix the renormalized param- 
eters, which are the physically relevant ones, then the corresponding bare 
parameters may be expressed relatively simply: 



n 2 



7o 




(12) 



From ( |T2| ) one can see that if T >> Q >> 7, then the same inequalities 
will hold for the bare quantities as well, and the proportional differences 
between the two sets of frequencies will all be of order (^) 2 , (p) 2 , or ^. The 
results in this subsection 2. 1 are all exact, and do not assume any particular 
relationships between the three frequencies. In the remainder of this paper, 
however, we will be interested in the case of extreme underdamping, with 
high cut-off frequency. We will therefore set 

7 = eVt 

n = eT , (13) 

where e is small. We will also consider -e(l+ln e) to be negligible. Only under 
this stronger assumption will the oscillator's final state be approximately 
pure0. 

The mode functions u^ and v u are both even functions in x. (Since the 
oscillator couples only to the even modes of the field, the odd mode functions 
are merely the usual sines. In our scenario, the odd modes will simply remain 
in their ground states forever, and so we will not refer to them explicitly 
hereafter.) For uj far from Qq, u^ and v u are essentially cosines; but they 
are distorted near the origin for u close to Qq, as one might expect. They 
possess several properties analogous to the orthonormality of cosines: 

dx Vu^Uu' (x) = tt5(u — lj') + q(u)p(uj') 



00 

00 7f 



du v u (x)u u (y) = -[5(x - y) + 5{x + y)] 
2 

du q(uj)uu(x) = / du p(uj)v UJ (x) = 

ii JO 

dui p{ui)q{ui) = —n . (14) 



These relations may all be verified by contour integration. The compu- 
tation is made easier if we use fll2"D to re-write ( |10|) as 



u u (x) 



Unj(X) 



dx F(x — x^Uuix') 



ftSrWr?-^- 2 



'QJ-O 



W 



([(fi 2 - uj 2 )(Yl + cu 2 ) + 2 7o r cj 2 ] coscux 

— 27 rQu;sinu;|a;|j (15) 



C(w) 



x 3J ([ w - (Q + ij)][u + (w - 17)] [w - ir][w + i(r + 2 7 )]e^ |x| 



Using (14) we can invert the transformation (g), to obtain 

1 



A 



n 



dx 4>(x)v UJ (x) — q(ui)P 
dx7r(x)u u (x) — p(uj)Q 



(16) 



It is then straightforward to verify from the standard commutation rela- 
tions of <p and 7T, and Q and P, that (jiff ) implies the canonical relations 



K,A^'] = 
[n£n£] = 

K,fl^] = i<J(u;-u/) 



Furthermore, we have 



A 1 /"OO . . 

ff=-y <z w ((itf) 2 + ^ 2 a 2 +) + h oM 



(17) 



(18) 



where Podd contains the unimportant odd mode operators B w and ITj . This 
demonstrates that @ is indeed equivalent to a free massless field. It is of 
course obvious that any quadratic model is equivalent to some spectrum of 
decoupled oscillators; but the fact that the single oscillator simply disappears 
like a drop in the continuous bucket of field modes, and does not alter the 
spectral density at all, is not so trivial. 



9 



2.2 Initial state 

Previous investigations of decoherence have typically considered direct prod- 
uct initial states, of the form 

|*i) = \0) field\lpa) oscillator , (19) 

where \Q) field is the vacuum state of the free field, and \ip a ) is a bimodal 
oscillator state, such as 

|Va> = c + e-^|0) osc + c_e^|0U , (20) 

|0)osc being the free oscillator ground state. For calculational convenience, 
we will instead use the initial state 

|*i> = c+e- iap \0) + c-e tap \Q) , (21) 

where |0) is the true, interacting ground state of the field-oscillator system. 
Since our objective in this paper is to study cases in which the initial 
oscillator entropy is effectively zero, we must show that the field-oscillator 



entanglement in the initial state (21 ) leads only to negligible initial entropy 
when the field is traced out. Since we have diagonalized the full Hamiltonian, 
the wave functional of the interacting ground state |0) in the rpj variables is 
simply a product of Gaussians: 

(n^Hz-vU 00 ^^) 2 . (22) 

Z is a normalization constant; a state (n^|0) is an eigenstate of the 11^ 
operators defined in (|16|)pl 



We can readily obtain from ( |2"2"D the corresponding reduced density matrix 
for the oscillator, 

Po(Q, Q') = JVn (Q, tt|0) (0\Q', vr) , (23) 

where we write T>ir for the measure to indicate that we want the limit where 
the integral is continuously infinite dimensional. Since according to (|1^) the 



"Both entities are of course only well-defined if we consider the continuous and infi- 
nite spectrum of oscillators to be the limit of a sequence of systems of finite numbers of 
oscillators. 

10 



operators IT^ are linear combinations of the operators n(x) and Q, the states 
|7r)|Q) are just different labellings of the states |II^). We therefore already 
have in equation (|22p the integrand of (P^D, and we need only now express 
the path-like integral over n in terms of the 11^ variables. 

We can deduce a straightforward way to do this from Equation (fTBl) . We 
first change variables in fl23|) by substituting 



. ir_, . -.1 P&) 



ir>(x), Q] - n^[7r(x)] - ^ Q , (24) 

'7T 



in (E2T), where 



- 4 1 f°° 

n^ = — p= / dx«a,(a;)7r(a;) . (25) 



Equation (|14T) then implies 



/ da;g(u;)n (J = 0. (26) 

Jo 



To integrate over tc(x), then, we will integrate over IL.,, with a delta-function 
inserted to enforce ([26|) and thus remove the Q sector from the integral: 

(n^-^Q) 2 +[(n^-^Q') 2 ' 



^ VJ TLV" W ^ 



Po(Q,Q') = Jd\ JvTle^^^^^e Uo °° 

= N e ~h MQ (Q 2+ Q' 2 ) e - AMn (Q~Q') 2 . (27) 

Here we introduce the renormalized mass M and ground state entanglement 
parameter A, which with our choice of frequency ratios ([13|) are given by 

M = M x (l + -J + 0(e 2 ' 

V 7T \l 

A = -l e (l + l ne ) + 0((e 2 ). (28) 

Since we have assumed (in order to have a pure state at late times) that A 
is negligibly small, we now observe that (|27D is effectively equal, for all Q,Q', 
to the density matrix of the (pure) ground state of an oscillator with mass M 
and natural frequency Q. Since the translation operators e ±mP in (pID only 
shift Q and Q' in fl2"T|) by ±a, it is also true that the reduced density matrix 

11 



formed from ([H]) differs negligibly from the one derived from ([19]). We have 
therefore shown that our oscillator has negligible entropy, when the field is 
traced out of our chosen initial state. In fact, our choice of (^TJ) instead of 
(|19"D will have no significant effect on our results (because we assume such 
weak coupling). 

For the remainder of this paper, we will set MQ = 1, and assume that, 
in the choice of units this implies, a 2 is large (of order e _1 ). This will mean 
that, even though the oscillator is very weakly coupled to the field, and has 
a dissipation timescale much longer than the dynamical timescale, the two 
Gaussians that are superposed in the initial state are far enough apart from 
each other that significant decoherence will occur on a very short timescale. 

2.3 Oscillator evolution 

By using the transition matrix in the momentum representation for a har- 
monic oscillator of frequency u, we easily obtain the wave functional of the 
final state into which the initial state (|2Tf) evolves at time t: 



(U A ,Ii B \^ t ) = iV 2 e 2jo " l 

x c + e^Jo <•> FK ' + c_e ^ Jo u y ' 1(29) 



We then use a technique like that employed in deriving (|27|) above to obtain 
the reduced density matrix. 

There will be four components: 

p(Q,Q';t) = |c + | 2 p++ + \c-\ 2 P— +c+c*_p+- +c^c\p_ + . (30) 

We find that the 'diagonal' components are given by 

P±±(Q,Q';t) = iVe-H^-mp+iQ'T-mp) 

xe Tias(t)(Q-Q>) 

xe- A '«-«'] 2 ; (31) 

while the cross-terms are 

p ±T (Q,Q';t) = jVe-Kl^^WI'+W'^WI 2 ) 

xe Tiaz(t)(Q+Q') 

xe -A[Q-Q' T 2y(t)] 2 

xe -a 2 ((l+4A)[l-^(t)]-z 2 (t)) _ ^ 32 ^ 

12 



Since A is negligible, we can ignore the last line of ( |31| ) and the second-last 
line of (J32|) . This leaves each of the four terms in p(Q, Q 1 ; t) as a separated 
function of the form 

p±±> = ^±±>(QW±AQ')- (33) 

The functions r(t) and s(t) are simple enough: 

r{t) = e~^[cosut-esmujt} + 0(e 2 ) 

s(t) = e-^[(l--e)smnt + 2ecosnt]-2ee- rt + 0(e 2 ). (34) 

IT 

It is clear from (|31|) and (|34|) that the diagonal terms describe Gaussian wave 
packets performing weakly damped oscillations. 

The functions y(t) and z(t), on the other hand, include some exponential- 
integral terms: 

y {t) = e-~ ft [{l--e-4A)cosQt-2esmQt]--e + 0{e 2 ) 

IT IT 

--e$t ((1 + i^)e iUt [Ei(-iQt) + id 

it \ 2 

+-e (e r4 Ei(-rt) + e~ r4 Ei(rt)) 
z (t) = e" 7 * [sin Qt + e cos Qt] + 0{e 2 ) 

--e3 ((1 + ittt)e int [Ei(-iQt) + m}) . (35) 

(The exponential-integral functions of imaginary argument appear together 
with iiT because we actually need a different branch of the Ei function than 
the standard one.) 

The cross-terms are very rapidly suppressed, because the exponent of the 
last term in (|32| 



D(t) = a 2 (l + 4A)[1 - y 2 (t)} - a 2 z 2 (t) , (36) 

grows on the cut-off timescale (see Figure 1), and a 2 is large. We emphasize 
that this decoherence occurs even in the extremely underdamped limit where 
A — > 0, and even when the initial state is not an exact direct product. 

The decohering factor D(t) grows rapidly because at times much less than 
a dynamical time f2 _1 , y(t) is approximately given by 

y{t) ~ 1 + -e (e ri Ei(-H) + e~ rt Ei(n) - 2C^ er - 2 ln(R)) + 0(e 2 ) . (37) 



IT 
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This function drops from unity on the cut-off timescale T" 1 . We can there- 
fore see that y(t) differs from r(t), and z(t) differs from s(t), on the same 
timescale as D(t) suppresses the crossterms. Hence there is a very short, 
very early time interval during which the four wave functions ip±±' appear- 
ing in p(Q, Q'; t) are all significant and distinct. Orthogonality, however, as 
opposed to mere distinctness, is what will be important for determining the 
eigenvalues and hence the entropy of the reduced density matrix. By the 
time y(t) has diverged from r(t) enough that ip± T are effectively orthogonal 
to ip±± 1 D{t) is already large, and ip± T may be ignored. We will therefore be 
able to neglect the distinction between if)±± and ip± T , and approximate the 
density matrix by the simpler, two-state form 

p(Q,Q';t) = \c + \ 2 MQW + (Q') + \c-?^-{QW-{Q') 
+ e- D ^c + c*_MQW-(Q') 

+ e-^c_ C ;f(Q)^((?') ) (38) 

where 

ip±{Q) = i>±±{Q) = Nh-^ QTr(t)]2Tias(t)Q . (39) 

2.4 Entropy evolution 

We can explicitly diagonalize p(t) in the approximation that ( |38D is valid by 
assuming that the eigenvector wave functions are of the form 

A (Q) = 5>±(A)^±(Q) . (40) 

± 

We then solve for the co-efficients A± by requiring them to be elements of 
the eigenvectors of the matrix 

/ \c+\ 2 + c+c ^ e -^W+« 2 (- 2 +- 2 ) | c+ |2 e -^W + c + c*_e- a ^ r2+s ^ \ 

1 |c_| 2 e~ D(4) +C_C* e~ a2(r2+s2) |C_| 2 + C-C* e ~D(t)+a 2 (r 2 +s 2 ) J ■ 

(41) 
The quadratic characteristic equation for M. yields the two eigenvalues 



^-\ 



1 ± a/1 - 4|c+| 2 |c_| 2 (l - e - 2a2e " 27t - e~ 2D ^] 



+ 0(e 2 ) , (42) 



when we take advantage of the facts that |c+| 2 + |c_| 2 = 1, e a ( r + s ) 

e _ a 2 e -2 7t &nd e -D(t) e -a 2 e" 
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We now specialize to the case c+ = c_ = 4=, where we have 



n/2' 

A ± ^i(l±[e- 2 ^ + e-^)]), (43) 

since e~ a e 7 e~ D ^> is extremely small for any t. 

The entropy S(t) for this case is plotted in Figure 2. It initially rises on the 
decoherence timescale from its initial negligible value to In 2, where it persists 
for many dynamical times, before declining on the dissipative timescale. The 
re-establishment of the purity of the oscillator state is clearly due to the fact 
that after a time on the order of 7 _1 In a the two shifted gaussians have lost 
so much amplitude that they begin to overlap and become indistinguishable. 
It then becomes less and less true that the oscillator is in a mixture of two 
orthogonal states. In the limit of complete relaxation, the ground state is 
reached, and this is of course a pure state. 

We can assess the accuracy of the approximation fl3~8|) by finding the 
eigenvalues of p{t) using ( |33"D instead. In this case we would have a fourth 
rank matrix in the analogue of (f|l]), and we would find four eigenvalues. Only 
two of these would be non- negligible, however, and they would turn out to 
differ insignificantly from (f4~2|). The eigenvectors we would find by this more 
accurate technique would differ somewhat, at very early times, from the two 
4>\ implied by ( |40| ) and ([11]). At these early times, the more exact eigenvectors 
would also include some non- vanishing amplitudes for the states represented 
by ?p± T (Q), with their dependencies on y(t) and z{t) instead of r(t) and s(t). 

2.5 Field evolution 

We now wish to consider the massless scalar field as the observed system, 
and to trace out the harmonic oscillator as unobserved. Since the initial 
state of the total system is pure, the non-zero eigenvalues of the reduced 
density matrix of the field are the same as those of the reduced density 
matrix of the oscillator. The entropy of the field is therefore the same as the 
oscillator entropy discussed in the preceding subsection. The problem that 
still remains, and which did not arise for the oscillator, is that of assessing 
where in spacetime the information associated with this entropy may be said 
to reside. 

There are of course many possible definitions of the term "information" , 
but for the purposes of this paper we shall consider that the information 
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problem will be solved by identifying the field state as a simple mixture of 
states that can be created from the vacuum by external sources. The external 
sources will be functions in spacetime, and the required information will be 
considered to reside in the regions where these sources have support. 

We will present the reduced density matrix R(t) for the field in the basis 
of field operator eigenstates, and so we will need to transform the field- 
oscillator state (ITj, II^ \ty t ) of Equation fl29|) from the IEJ, IT^ representation 
into the A w , B^ representation. We then invoke one of the eigenvalue relations 
corresponding to Equation ( |16|) 



A 



<ix0(x)f aJ (x) — q{ui)P 



and trace over P to obtain the reduced density matrix 

R[<j>, 4>'; t] = JdP (A u [<f>, P]\%) (%\A„[<f>>, P)) 



(44) 



(45) 



The reduced density matrix for the field will again be the sum of four 
contributions, 



where 



R[<f>,<f/;t] 

R±± = 

R± T = 



\c+\ 2 R ++ + |c l 2 i?__ + c + c*_R, + c„cli?__ 



(46) 



*±[<t>;t]n[<t>';t}e 

* ± [<f>;tm[<f>';t}e 



f°° dxU+<f>')L(x) 

•J — oo 

P° dx((j,+<t>')L{x) 



Here we have defined 
L(x) = 



xe -a\r 2 +s 2 ) e ±a[ir(t)~s(t)]f X oo dx^+cf>>)L(x) 



1 f°° 

— / duo uj q((jj)vu{x) 
n Jo 



(47) 



* \o(x)-t] = z-* e -ur d "K\^ sin " t+aq ( w >^-j-oc dx(t,( - x >"^ 



(x);t\ 



±ia I dx<b(x)K(x,t) 

X e J -°° 
The function K(x, t) will be defined below. 



(48) 
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As in (|38|), we can simplify fl47j ) by approximating it as 

R±± ^ *±[0;t]*i[0';t] 

i? ±T ~ f ± [^;t]^^t]xe- aV211 . (49) 

This approximation can evidently lead to large errors in evaluating expecta- 
tion values of 4>(x) for x close to the origin, where L(x) is non- negligible; but 



as with the analogous approximation in (|38|) , it retains the actual behavior 



of the entropy very well. Equation ([EJ) and the inner product 

<*_(W+W> = e- D(4) (50) 



imply that the eigenvalues of the density matrix (|46|) are again given by 
as should be the case. 

In (^), there are only two wave functionals \I/± that characterize the 
state of the field. The states they represent may be created by an external 
source, linearly coupled to both </> and ir: 

where the time-independent state |\&o) is simply |\l/±) with a set equal to zero. 
(This state is not precisely the vacuum state of the field, but it resembles 
it closely except near x = 0, where it has been polarized by the unobserved 
oscillator.) 

J and K are the external sources which describe the spacetime location 
of the information: 

J(x,t) = — dx' F(x — x')r(t — \x'\) 
2 J-t 

K(x,t) = -(2r(t)F(x)-F(x + t)-F(x-t) 

-(1 + -e) f dx' F(x - x')s(t - \x'\)\ . (52) 

From ([52]) it is clear that in the late time limit, all the information has 
propagated away from the origin. In this late time limit e~ a e — > 1, and 
the state of the field, with the oscillator traced out, becomes pure. 

Note that the behavior of J and K is not particularly sensitive to the 
small amounts of initial field-oscillator entanglement that are at issue in 
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choosing (|21| ) instead of a direct product state. The operator exponent in 
(|5TD is simply the projection, onto the field sector of the total Hilbert space, 
of the time-evolved P operator: 

P{t) = r{t)P - Mtts{t)Q + / dx [J{x, t)vr(x) + K{x, t)<j>{x)\ . (53) 

If the information initially in the oscillator is characterized by P, then the 
information propagation into the field will be described (up to ambiguities 
near the origin) by fl52|). This clearly shows that the information which was 
characterized by the initial value of P of the oscillator will become non-local 
in the late time limit. 



3 Discussion 

The evolution of the oscillator-field system shows that the entanglement en- 
tropy is roughly constant until the oscillator approaches its ground state. 
Moreover we have shown that the information is hidden in a very non-local 
way. 

Both J(x, t) and K(x, t) are important in characterizing the state of the 
field. J(x, t) behaves very much like classical radiation from a source, and 
implies that some of the information initially in the oscillator propagates into 
the field in the same way that one might naively expect. In contrast, K(x, t) 
possesses sharp spikes, whose width is on the cut-off scale, which propagate 
away or decay in place. The propagating spikes are evidently the couriers 
for the rapid shedding of information that is associated with decoherence. 
The spike in K(x, t) which remains at x = decays only on the dissipative 
timescale. This could be interpreted as describing information which remains 
in the 'quantum hair' of the oscillator, and is only slowly radiated away. It 
should be pointed out, however, that our approximation in (f49|) breaks down 
near the origin. 

Furthermore, assigning the degrees of freedom near the origin to the field 
or to the unobserved oscillator is to a large degree arbitrary, since unitary 
transformations that alter such assignments can leave the rest of the model 
essentially unchanged. It could well be that the information pool described 
by the non-propagating term in K(x, t) is best considered as belonging to 
a 'dressed' version of the oscillator. This type of ambiguity seems likely to 



be a common feature of problems in which information loss and ultraviolet 
regulators are closely related, since one does not expect a uniquely specified 
ultraviolet cut-off to exist. In fact, the best procedure in such cases might 
be to search for an optimal kind of UV regulation based on the location of 
information. 

We can arrange to present the state of the field (with the oscillator treated 
as unobserved) in the form of a Wigner function. We can do this by not- 
ing that (i) the field state is effectively a mixture of pure states |\l/±(£)), 
and (ii) at any fixed time t, the states |\I/±(£)) belong to a subspace of the 
field's Hilbert space which may be mapped onto the Hilbert space of a single 
harmonic oscillator. The second assertion is justified by the inner product 
Equation fl5H|), which is the correct inner product for two coherent states with 

annihilation operator eigenvalues a± = ±y — ^.[] Therefore, by applying to 
and (H3) the mapping 



|* ± (t)> - | ± sjD(t)/2) osc , (54) 

the state of the field may be described at any given time by a mixture of 
single-oscillator coherent states \a±) osc . 

Wigner functions for such mixed states are easily calculated. The function 
for the field at t = is single gaussian corresponding to the ground state. 
The Wigner functions which represent the state of the field at several later 
instants are plotted in Figures 3.1 - 3.3. Figure 3.1 shows the state of the 
field very soon after t — 0, as it is just getting excited into a mixture of two 
coherent states. The two overlapping gaussian peaks that can be discerned in 
3.1 will separate on the decoherence timescale. (No rapid growth in energy 
is associated with this sudden separation: it is only the inner products of 
the field states, and not their energies, that fit the analogy with oscillator 
coherent states.) Figure 3.2 describes the field during the long interemediate 
epoch of its evolution: a mixture of two well-separated gaussians. Finally, 
when the oscillator ends up in its ground state, the field regains its purity. 
The rapid oscillations near the origin in Figure 3.3 are symptoms of quantum 
coherence. 



3 By repeating the analysis of Section 2 starting with more general initial states, it 
is straightforward to define a class of field states which represent, in the sense we are 
considering, the full complex plane of coherent states, with complex eigenvalues a. 
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These figures clearly show the loss and eventual restoration of quantum 
coherence in the field. In order to determine the degree of purity of a har- 
monic oscillator, however, one must have some means of measuring at least 
the gross features of its Wigner function. In our case, the observables which 
must thus be measured (in whatever combination) are the non-local oper- 
ators J^° 00 dx[J(x,t)7r(x) + K(x,t)(j)(x)} and its canonical conjugate. Since 
these operators are non-local on the dissipation scale, it would seem to be 
very difficult to actually observe the asymptotic purity of the state of the 
field. 



This system can be thought as a specific example of Page'spO[] alternative 
outcome of black hole evaporation. He suggests that information might get 
out of black holes through radiation, and has shown that the information 
might not show up in an analysis perturbative in M P i anck /M. Our system 
seems to behave in this way, as in our case the information is not recovered 
until the very end of the decay of the oscillator. If we were to start the oscil- 
lator in a mixture of the two gaussian states, instead of in the superposition 
we have discussed, we would not notice any difference by examining the field 
until the oscillator had relaxed to very near its ground state. How long this 
takes to occur, in our model, is dependent on the initial state. 

We should also point out that we do not find a strict relation between 
the energy of the oscillator and the rate of information exchange (one of the 
folkloric statements used as an argument for the loss of coherence in black 
hole radiation). Thus it might well be possible that black hole evolution does 
preserve quantum coherence, but that it is very hard, if not impossible, to 
recover all the initial information in a set of local observations of the final 
state. 
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Figure captions. 

Figure 1. Behavior of D(t) of equation (36), with the parameter a set equal 

to 10. D(t) is the term responsible for the suppression of the interference 
term of the two-gaussian initial state; it increases rapidly (on the decoherence 
timescale) and remains large. 



Figure 2. Evolution of the entropy of the oscillator for three initial states of 

different values of the initial separation parameter a. The entropy increases 
rapidly during the decoherence phase and remain esentially constant until 
the oscillator enters the final stages of its decay towards the ground state. 
Since this occurs only when the oscillator has lost all but a fixed amount of 
its initial energy, the length of time before recoherence depends on the initial 
state. 



Figure 3. Evolution of the the degree of freedom of the field which is excited 

by the oscillator (i.e., the mode spanning linear combinations of the time- 
dependent, non-locally excited states \^±(t))). It starts in its ground state, 
is excited into a mixed state by the oscillator, and approaches a pure excited 
state as the oscillator decays to its own ground state. The sharp oscillations 
in 3.3 are witnesses to the purity of the field state at late times. 
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This figure "fig2-l.png" is available in "png" format from: 
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